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BACKWATER FUNCTIONS BY NUMERICAL 
INTEGRATION 


By CLINT J. KEIFER,' A. M. ASCE, AND 
HENRY HSIEN CHU* 


1 Asst. Engr. of Sewer Design, Dept. of Publie Works, Bureau of Eng., Chicago, UL. 
2 Senior Asst. Sewer Designer, Dept. of Publie Works, Bureau of Eng., Chicage, Ill. 


Synopsis 
Because of the importance of nonuniform flow and the particular difficulties 
in computing backwater curves in closed conduits, the authors have devised 
new procedures facilitating backwater computations. Complete derivations 
of equations for computing backwater functions are presented, and adaptations 
of the step method are suggested as additional aids in this field. 


INTRODUCTION 


Notation.—The letter symbols introduced in this paper are defined where 
they first appear, in the text or by illustration, and are assembled alphabetically 
in the Appendix for convenience of reference. 

Hydraulics engineers are repeatedly required to compute backwater curves. 


In a natural water course, the shape of the cross section may vary greatly 
among sections. This necessitates use of one of the many step methods. If 
the channel is uniform, the computations may be minimized through the 
integration method appropriate for the particular problem. For channels that 
are wide in comparison to the depth of flow, for example, one may choose the 
Bresse backwater functions. The method most applicable to various shapes 
of open channel was devised by Boris A. Bakhmeteff,? Hon. M. ASCE, and 


- . “Hydraulics of Open Channels,” by Boris A. Bakhmeteff, McGraw-Hill Book Co., Inc., New York, 
¥., 2088. 


supplemented to in¢lude velocity head by M. E. Von Seggern,‘ A.M. ASCE. 


‘Integrating the Equation of Nonuniform Flow," by M. E. Von Seggern, Transactions, ASCE, Vol. 
115, 1950, p. 71. 


These methods cannot be applied accurately to closed conduits when flow 
is near the top. The writers propose a method for separating many of the in- 
volved factors in the nonuniform flow equation so’ that it can be integrated 
numerically with respect to the relative depth of flow, y/D, in which y represents 
the depth of flow and D is the diameter of the conduit. Circular sections have 
been taken as examples, and tables corresponding to the Bresse function for 
shallow rectangular sections have been prepared. The application of this 
circular section is limited, but the same principle and procedure can be applied 
to any type of conduit or channel. 


t : 
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NUMERICAL INTEGRATION 


Equation for Nonuniform Steady Flow.—The equation for nonuniform steady 
flow can be written in the following form: 


in which S is the energy gradient, x represents distance (measured upstream), 
V is the average velocity at a section, g is the acceleration produced by gravity, 
and S, denotes the slope of the invert or bed. This equation indicates that the 
slope of the energy gradient is equal to the change in the depth of flow plus the 
velocity head, with respect to a unit length (dr) plus the bottom slope. 
Substituting Q?/A? for V? (Q representing discharge and A being the cross- 


sectional area of the flow) S = dy d ( Q ) + So. By differentiation, 


dz * dx\2g A? 


Eq. 3 is the general differential equation, taking So positive when it is in- 
clined downstream and z is measured upstream. Because the expression be- 
tween the parentheses in Eq. 3 is a function of depth, this equation may be 
written: 

dx = f(y) dy 


In this form it may be integrated graphically if the size and shape of the 
conduit or channel, the invert slope, and the discharge are held constant. 
Each time one of the parameters is changed, the graphical integration must be 
performed again. 

Eq. 3 has not been completely evaluated mathematically without approxi- 
mations. Mr. Bakhmeteff has integrated Eq. 3 by assuming that the square 
of the conveyance at any depth of a channel varies as the ith power of the 
depth; so that, if k is the conveyance and 1 any positive number, k? « y'. He 
also neglected the effect of the change in velocity head. This procedure ne- 
cessitates the use of short steps so that the change in velocity head within each 
step would not be too great. Mr. Von Seggern makes the same assumptions 
as Mr. Bakhmeteff and introduces another approximation, A*/b, « y™, which 
takes account of the velocity-head changes, m being any positive number. 
These assumptions enabled both writers to integrate Eq. 3 analytically and to 
obtain approximate solutions for most types of open channels. The closed 
conduit presents a different situation in which the values of i and m, instead of 
being constant, vary rapidly when the water depth approaches the crown— 
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_dy_ @ dAdy 
dA 
Since the top width of water surface b,. = dy’ 
| 
{ 
| 


thus rendering the basic assumption inapplicable. Fig. 1 shows curves de- 
seribing this variation. The subscript F indicates a full section. The ex- 
pressions for the exponents t and m are 


a 


t= 2cot@ =2 


in which 6 and @ are slopes of the curves in Fig. 1. Values thus determined 
appear in Table 1. Four parameters are present in the integrand of Eq. 4. A 


TABLE 1.—VaRIaTION OF t AND m IN CIRCULAR CONDUITS 


0.30 | 040 0.50 0.60 0.65 0.70 0.75 0.80 0.85 0.90 
3.12 3.01 2.80 2.47 1.93 1.20 0.85 
3.87 3.87 3.87 3.89 3.99 4.01 4.14 4.42 4.92 5.98 


method has been devised for removing two of the parameters from under the 
integral sign so that for a specified conduit shape the numerical integration 
must be performed only over the range of one remaining parameter. The fol- 
lowing derivation thus required no approximations. 

If a particular shape of conduit is chosen—such as a circular se¢tion—and 
the dimensionless ratios y/D and x/D are used to express depth of flow and 
length along the channel, then Eq. 3 becomes 


The flow Qeap = kev So 
and 
Q = kVS = keV Sp = 
in which Q.ay represents the capacity of the conduit flowing full and having the 
energy gradient equal to the bottom slope, the subscript zero indicates the 
channel bottom, and the subscript o indicates normal depth. 


Therefore, 
(ir) =s (2) (#) 
A OF k Oo Qeap k 


Observing that (kr/k)? is a function only of y/D (designated f,(y/D)) and 
replacing Q/Qcap by a for convenience, it is concluded that 


a (Io 
ED 
and 
a (log £4* ) 
m = = ————_ 
(or ) 
—0.54 
- 
gf y 
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In Eq. 6, the value of o after multiplying by D»*, is also a function of 


the proportional depth y/D—written as 


in which the subscript ¢ corresponds to conditions at critical depth. The 
graphs of f, (y/D) and fz (y/D) for circular conduits are shown in Fig. 2. Sub- 
stituting Eqs. 9 and 10a in Eq. 6, 


a’? 


So Dt » 
(5) 


in which y, and y. represent flow depths at two consecutive specified sections. 
Iq. 116 can be integrated by assigning values to a and integrating numerical 


1 
values of —————————- and | at each increment of 


0.01 in the value of (y/D), using Simpson’s rule. The integration can be 
started at any reference value of (y/D), such as zero or any convenient depth. 
Designating members of Eq. 116 by 


Also, 
| 
and | 
- 
| 
| 
and 
y ) - 
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Eq. 116 becomes 


D Q? 
For a certain shape of section, the values of ® and ¥ are functions only of 
the proportional depth anda. They shall be called backwater functions. These 
functions have been computed for circular conduits and are exhibited in both 
graphical and tabular form in Figs. 3 and 4 and in Tables 2 and 3. 


‘ More extensive and detailed tables of the backwater functions may be obtained free upon request, 
from the authors. Address requests to Messrs. Keifer and Chu, Department of Publie Works, 188 West 


Randolph Street, Chicago 6, Il. 


TABLE 2.—BackwatTerR Function ® For CircuLar SECTIONS 


VALUES OF a 


0.90 1.00 ‘ 1.20 1.40 1.60 
0 0.9166 | 0.4075 | 0.2543 
0.8567 | 0.3830 | 0.2397 
0.7822 | 0.3555 | 0.2238 
0.7016 | 0.3269 | 0.2074 
0.6203 | 0.2982 | 0.1909 
2.7491 | 0.5427 | 0.2701 | 0.1748 
1.4754 | 0.4714 | 0.2432 | 0.1591 
1.0082 | 0.4076 | 0.2178 | 0.1441 
0.7486 | 0.3515 | 0.1942 | 0.1299 
0.5803 | 0.3025 | 0.1724 | 0.1166 
0.8911 | 0.4619 | 0.2600 | 0.1524 | 0.1043 
0.6046 | 0.3740 | 0.2233 | 0.1343 | 0.0928 
0.4518 | 0.3065 | 0.1916 | 0.1179 | 0.0823 
1.5967 | 0.3521 | 0.2533 | 0.1643 | 0.1031 | 0.0727 
0.5844 | 0.2810 | 0.2105 | 0.1407 | 0.0899 | 0.0639 
0.4016 | 0.2275 | 0.1765 | 0.1202 | 0.0781 | 0.0560 
0.3001 | 0.1858 | 0.1465 | 0.1026 | 0.0676 | 0.0488 
0.5102 | 0.2329 | 0.1527 | 0.1228 | 0.0873 | 0.0583 | 0.0424 
0.3311 | 0.1845 | 0.1259 | 0.1023 | 0.0741 | 0.0500 | 0.0366 
0.2403 | 0.1479 | 0.1039 | 0.0857 | 0.0627 | 0.0428 | 0.0315 
0.4083 | 0.1825 | 0.1196 | 0.0858 | 0.0713 | 0.0529 | 0.0364 | 0.0269 
0.2596 | 0.1416 | 0.0970 | 0.0708 | 0.0593 | 0.0444 | 0.0308 | 0.0229 
0.1851 | 0.1113 | 0.0789 | 0.0583 | 0.0492 | 0.0371 | 0.0259 | 0.0193 
0.1380 | 0.0882 | 0.0640 | 0.0478 | 0.0406 | 0.0309 | 0.0217 | 0.0163 
0.1971 | 0.1051 | 0.0701 | 0.0512 | 0.0391 | 0.0334 | 0.0255 | 0.0180 | 0.0136 
0.1385 | 0.0811 | 0.0557 | 0.0410 | 0.0318 | 0.0274 | 0.0210 | 0.0149 | 0.0112 
0.1014 | 0.0628 | 0.0442 | 0.0331 | 0.0257 | 0.0221 | 0.0171 | 0.0122 | 0.0092 
0.0756 | 0.0487 | 0.0347 | 0.0265 | 0.0206 | 0.0178 | 0.0138 | 0.0099 | 9.0075 
0.0569 | 0.0377 | 0.0273 | 0.0210 | 0.0164 | 0.0142 | 0.0110 | 0.0080 | 0.0061 
0.0430 | 0.0291 | 0.0212 | 0.0164 | 0.0130 | 0.0113 | 0.0088 | 0.0063 | 0.0048 
0.0324 | 0.0223 | 0.0165 | 0.0127 | 6.0101 | 0.0088 | 0.0069 | 0.0050 | 0.0038 
0.0155 | 0.0110 | 0.0082 | 0.0064 | 0.0052 | 0.0045 | 0.0035 | 0.0026 | 0.0020 
0.0068 | 0.0049 | 0.0038 | 6.0029 | 0.0024 | 0.0021 | 0.0016 | 0.0012 | 0.0009 
0.0027 | 0.0019 | 0.0015 | 0.0012 | 0.0009 | 0.0009 | 0.0006 | 0.0005 | 0.0004 
0.0008 | 0.0006 | 0.0005. | 0.0004 | 0.0003 | 0.0003 | 0.0002 | 0.0002 | 0.0001 

0.0002 | 0.0001 | 0.0001 | 0.0001 | 0.0001 | 0.0001 _ —_ — 
0 0 0 0 0 0 0 


| 
| 
| Ratio, 
D 0.40 
1.00 | 0 
0.98 | 0.0235 
0.96 | 0.0468 
0.94 | 0.0700 
0.92 | 0.0932 
0.90 | 0.1165 
0.88 | 0.1398 
0.86 | 0.1632 
0.84 | 0.1868 
0.82 | 0.2105 
0.80 | 0.2344 
0.78 | 0.2585 
0.76 | 0.2830 
0.74 | 0.3077 
0.72 | 0.3329 
0.70 | 0.3585 
0.68 | 0.3848 
0.66 | 0.4117 
0.64 | 0.4304 
0.62 | 0.4683 ba 
0.60 | 0.4985 
0.58 | 0.5304 
0.56 | 0.5647 
0.54 | 0.6022 
0.52 | 0.6446 
0.50 | 0.6947 
0.48 | 0.7591 
0.46 | 0.8585 
0.44 | 3.3677 
0.42 | 0.1837 
0.40 | 0.1107 
0.35 | 0.0419 
0.30 | 0.0169 
0.25 | 0.0062 
0.20 | 0.0019 
0.15 | 0.0004 
0 0 
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TABLE 3.—Backwarer Funcrion For CIRCULAR SECTIONS 


VALUES OF a 


D 0.40 0.60 0.70 0.80 0.90 1.00 1.07 1.20 1.40 1.60 


0.0245 | 0.1128 | 0.0680 | 0.0482 


1.00 0 0 0 0 0 
0.98 0.0003 | 0.0003 | 0.0004 | 0.0006 | 0.0009 | 0.0016 | 0.0186 | 0.1120 | 0.0677 | 0.0479 


0.96 0.0008 | 0.0010 | 0.0012 | 0.0015 | 0.0024 | 0.0053 0 0.1105 | 0.0671 | 0.0477 : 
0.94 0.0015 | 0.0019 | 0.0023 | 0.0028 | 0.0043 | 0.0098 | 0.0824 F 0.1079 | 0.0663 | 0.0470 | 


0.92 0.0023 | 0.0029 | 0.0036 | 0.0045 | 0.0067 | 0.0150 | 0.1614 JF 0.1051 | 0.6052 | 0.0466 
0.90 | 0.0034 | 0.0042 | 0.0050 | 0.0063 | 0.0096 | 0.0215 [0.2500 | 0.1028 | 0.0641 | 0.0458 
0.88 0.0044 | 0.0055 | 0.0067 | 0.0086 | 0.0132 | 0.0303 | 0.1852 | 0.0983 | 0.0627 | 0.0452 
0.86 0.0058 | 0.0072 | 0.0087 | 0.0112 | 0.0175 | 0.0438 | 0.1604 | 0.0950 | 0.0613 | 0.0444 
0.84 0.0071 | 0.0090 | 0.0110 | 0.0143 | 0.0230 | 0.0667 | 0.1439 | 0.0913 | 0.0598 | 0.0435 
0.82 0.0087 | 0.0112 | 0.0136 | 0.0180 | 0.0299 0.1321 | 0.0879 | 0.0583 | 0.0426 
0.80 0.0106 | 0.0136 | 0.0166 | 0.0225 | 0.0389 | 0.1806 | 0.1230 | 0.0847 | 0.0568 | 0.0416 
0.78 0.0128 | 0.0164 | 0.0204 | 0.0282 | 0.0524 7 0.1546 | 0.1153 | 0.0816 | 0.0552 | 0.0405 
0.76 0.0150 | 0.0196 | 0.0247 | 0.0349 | 0.0768 | 0.1394 | 0.1088 | 0.0784 | 0.0537 | 0.0396 
0.74 0.0177 | 0.0236 | 0.0300 | 0.0440 7 1.8457 | 0.1288 | 0.1032 | 0.0755 | 0.0521 | 0.0384 
0.72 0.0206 | 0.0278 | 0.0361 | 0.0571 | 0.1950 | 0.1200 | 0.0980 | 0.0727 | 0.0505 | 0.0375 
0.70 0.0240 | 0.0333 | 0.0443 | 0.0817 | 0.1650 | 0.1130 | 0.0933 | 0.0699 | 0.0489 | 0.0363 
0.68 | 0.0280 | 0.0395 | 0.0556 | 0.1747 ] 0.1465 | 0.1066 | 0.0891 | 0.0673 | 0.0473 | 0.0354 : 
0.66 0.0324 | 0.0475 | 0.0724 [0.2320 | 0.1352 | 0.1011 | 0.0851 | 0.0648 | 0.0458 | 0.0343 | 
0.64 0.0377 | 0.0576 | 0.0978 | 0.1882 | 0.1263 | 0.0960 | 0.0813 | 0.0622 | 0.0442 | 0.0332 
0.62 0.0438 | 0.0721 | 0.2500 | 0.1680 | 0.1189 | 0.0916 | 0.0777 | 0.0599 | 0.0426 | 0.0322 
0.60 | 0.0511 | 0.0915 [0.2780 | 0.1545 | 0.1120 | 0.0870 | 0.0743 | 0.0575 | 0.0411 | 0.0310 
0.58 0.0600 | 0.1290 | 0.2265 | 0.1437 | 0.1054 | 0.0830 | 0.0710 | 0.0550 | 0.0396 | 0.0299 | 
0.56 0.0705 | 0.3909 | 0.2006 | 0.1343 | 0.0997 | 0.0790 | 0.0678 | 0.0529 | 0.0381 | 0.0288 i 
0.54 0.0840 | 0.3080 | 0.1828 | 0.1260 | 0.0944 | 0.0753 | 0.0647 | 0.0506 | 0.0365 | 0.0278 i 


0.52 | 0.1014 | 0.2580 | 0.1686 | 0.1186 | 0.0895 | 0.0717 | 0.0618 | 0.0485 | 0.0351 | 0.0266 
0.50 | 0.1253 | 0.2293 | 0.1562 | 0.1117 | 0.0850 | 0.0682 | 0.0588 | 0.0463 | 0.0336 | 0.0256 
0.48 | 0.1609 | 0.2085 | 0.1452 | 0.1054 | 0.0803 | 0.0648 | 0.0561 | 0.0441 | 0.0321 | 0.0244 
0.46 | 0.2235 | 0.1912 | 0.1350 | 0.0993 | 0.0761 | 0.0616 | 0.0533 | 0.0420 | 0.0306 | 0.0233 
1.4143 | 0.1770 | 0.1256 | 0.0935 | 0.0718 | 0.0582 | 0.0505 | 0.0399 | 0.0291 | 0.0221 
0.1643 | 0.1170 | 0.0880 | 0.0680 | 0.0552 | 0.0478 | 0.0378 | 0.0276 | 0.0210 
0.1528 | 0.1090 | 0.0826 | 0.0641 | 0.0519 | 0.0451 | 0.0357 | 0.0261 | 0.0199 
0.1271 | 0.0917 | 0.0700 | 0.0546 | 0.0444 | 0.0386 | 0.0306 | 0.0224 | 0.0171 


0.1038 | 0.0755 | 0.0580 | 0.0454 | 0.0370 | 0.0321 | 0.0255 | 0.0188 | 0.0143 
0.0825 | 0.0607 | 0.0466 | 0.0366 | 0.0297 | 0.0258 | 0.0205 | 0.0151 | 0.0114 
0.0627 | 0.0462 | 0.0356 | 0.0281 | 0.0227 | 0.0197 | 0.0157 | 0.0115 | 0.0088 
0.0442 | 0.0327 | 0.0250 | 0.0197 | 0.0158 | 0.0138 | 0.0110 | 0.0081 | 0.0062 


0 0 0 0 


0 0 


0 0 


0 


Referring to Eqs. 7, when flow is at normal depth, 


Therefore, a® f; (yo/D) — 1 = 0, and the integrand in Eq. 11a becomes infinity. 
Infinite length is required to make any change in normal water depth, meaning 
that the backwater curve approaches the normal depth asymptotically. Eq. 
10b indicates that, when flow is at critical depth, the numerator in Eq. lla 
vanishes. This indicates the length of backwater curve to be zero for any 
change in water depth. 

The heavy lines in Tables 2 and 3 show the location of normal depth. The 
tabulated values can be interpolated for intermediate values of ® and V, but not 
across the heavy lines, because a backwater can never pass through the normal 
depth. The bottom, y/D = 0, cannot be used as a starting point of integration 
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Ratio, 
& 
0.30 | 0.2432 | 
0.25 | 0.1903 
0.20 0.1425 4 
0.15 | 0.1000 
0 0 | | 


in evaluating ® and W for depth above the normal depth. The zero value of 
# and V in Tables 2 and 3 shows where the integration has been begun in each 
case. 

Eq. 12 indicates that the backwater functions ® and ¥ will apply to all 
cases, regardless of the diameter and discharge values. The usage of these 
functions, however, must be restricted to conduits having bottom slopes in the 
direction of flow either steeper or milder than the critical slope. If the bottom 
is horizontal, Sy = 0,a@ = «, and Eq. 9 becomes meaningless. In the case with 
adverse slope, So is negative, and 


Qeap Kp So 


Therefore, a? = (Q/Qeu))? is also negative. The integration process must be 
repeated, using these expressions to determine values of ® or V for the adverse 
slope. 

In the case of an open channel, any convenient linear measure (for example, 
the bottom width) can be used instead of the full depth of closed conduits. The 
remaining expressions can be derived in the same manner as for closed conduits. 

The preparing of the tables for backwater functions for one type of section 
is laborious—but justified if the section will be used frequently and if many 
backwater curves must be computed. The step method may be preferred, 
however, for hydraulic investigation of conduits having other cross sections 
that will not be used so often. 


BacKWATER CurRvES IN UNniIFORM CHANNELS BY THE STEP METHOD 


The step method for computing backwater curves has certain advantages 
over the integration method. It provides a complete profile of the water sur- 
face, values of the velocity at each section, and many other wanted factors. 

In the step method, Eq. 3 is 


A (y + =) 
= 
Sav So 

in which S,, is the average value of the slope of the energy gradient at the two 
ends of a reach. For computing the energy slope and velocity head, many 
laborious steps are involved in commonly used forms of backwater computation. 
To save this labor, the writers have devised two new approaches which are 
supplementary te the integration method. 

In the first of these, from Eq. 8, 


Sr 1 


Because Sp can be computed from given values of the discharge Q, and diameter 
or size D, values of S for any given water depth can be determined from Eq. 16. 
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Also for given values of discharge und diameter the velocity head can be com- 
puted from the relationship, 
an (2 
29 


in which h, is the velocity head, and hy is the velocity head of a full section. 

Eqs. 16 and 17 are represented graphically in Fig. 5 for circular sections, Sp 
having been determined by the Manning formula. Table 4 facilitates the 
computations of slope full and velocity head flowing full, as expressed in the 
following equations: 


Q 2 
Sp =a ( n= 0.013 
Q 


2 
(35) 


in which a and 6 are quantities evaluated in Table 4. Fig. 5 provides a scale to 


TABLE 4.—Vatues or Factors 1n Eas. 18 


COEFFICIENTS CorFFICIENTS 
Diameter, D, Diameter, D, 


| 
| 
| 
| 
| 


0.02450 
0.01696 
0.01202 
0.008701 
0.0064 10 
0.004806 
0.003656 
0.002819 
0.002199 
0.001735 
0.001383 
0.001112 
0.0009023 
0.0007376 
0.0006079 


0.03636 


aid in computing the critical depth, a more complete version of which has been 
published elsewhere.’ Detailed computation procedure is shown in Table 5. 


§“Hydraulies of Steady Plow in Open Channels,” by Sherman M. Woodward and Chesley J. Posey, 
John Wiley & Sons, Ine., New York, N. Y., 149, p. 42. 


1 
in Feet 4 in Feet 
| a | a 
4.822 
3.837 
3.091 
2.517 
2.072 
1.719 j 
1.439 
1.214 
1.031 j 
0.8814 
0.7578 
0.6553 
383-8 


= put = ag ‘poyjeur Ut ost 40j QT “bey = PUB ‘oas sad NI = % uy» 


1Z0't 
eye's G88 
009 
690'L 

CHL 

68 


tic 


cel 

cel 

cel 

Chl OFl 

ost 9 “M4 
£1z Col j 


d 


t 4 3 260 000'T 
i £260 2160 
siv0 L660 299 0 4 AIS yf 
£100 coro 09 0 
c6e°0 LE 


aD (or) 


| | 
| 
| | 
| 
| peer 
| 
| 
| 
| 
| 
RON 
| 
! 
| 
| 
t 
| 
| 
S| 
| | 
| 
| 
=| 
=| 
| 
| 
| | 
| | 
| 
383-9 
| 


Another method used for obtaining the energy slope and velocity head is the 
preparation of charts such as Fig. 6 for cireular conduits. In Fig. 6, the heavy 
broken lines with arrowheads show the operations when the size, discharge, and 
depth of flow are specified. Backwater curves can be computed by the pro- 
cedure shown in the following example (Table 5). 


ILLUSTRATIVE EXAMPLE 


A circular conduit of 72-in. diameter has a discharge of 160 cu ft per sec, a 
bottom slope of 1 ft per 1,000 ft and-a Manning friction factor of n = 0.013. 
The distance is to be determined between the control section at critical depth 
and the section at which the conduit becomes full: The discharge capacity 


Qeap equals 1s Ay Rip S'y or 134 cu ft per see (with Rp denoting the 
hydraulie radius of ow ‘nanan Since Q/D! = 1.82, the critical depth is 
obtained from Fig. 5 as * 5 = 0.575, or y. = 3.45. The type of backwater 
curve is M, in which the water depth y is less than y, but greater than y,. 
Integration Method. —From Eq. 12, = (2 #1) — 


Also, a = = 1.2, 6,000 ft, and 19.960 Thus, 21.2: = 
5,245 — 1,166 = 4,079 ft. 

Step Methods—The computations using the first step method are shown in 
Table 5 with the aid of Fig. 5. The distance is «;.,2 = 4,038 ft. The computa- 
tions using the second step method appear in Table 5 with the aid of Fig. 6. 


This procedure yields a value of x;.2 equal to 4,021 ft. 


CONCLUSIONS 


1. Most natural or artificial channels contain nonuniform flow. Uniform 
flow, which is usually assumed in making hydraulic computations, is actually 
rare. New methods are presented for solving nonunifiorm flow in prismatic 
channels, and the almost universal presence of nonuniform flow is emphasized. 

2. Tables of backwater functions ® and W¥ for circular conduits are presented 
in this paper. The computation of backwater curves is considerably shortened 
by the use of these tables. as demonstrated by the example. Tables could be 
prepared as needed for any other types of conduits. 

3. The time-consuming step method computations of the energy gradi- 
ent and the velocity head can be simplified by the two devices suggested in this 
paper. The writers believe these to have been heretofore unpublished. Both 
procedures require the preparation of charts before hand. Fig. 5 will yield more 
accurate results than Fig. 6 although values in Fig. 5 are expressed as a percent- 
age of those of the entire section, and one step of conversion is needed to deter- 
mine the actual values. Furthermore, more labor is required in preparing Fig. 
6 than Fig. 5. Therefore, the first method seems preferable. Charts similar to 
those used with the first method can also be prepared for open-channel sections; 
any convenient linear measure may be substituted for depth full in closed con- 
duits (for example, bottom width). 


] 
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APPENDIX. NOTATION 


The following symbols, adopted for use in the paper, and for the guidance 
of discussers, conform essentially to American Standard Letter Symbols for 
Hydraulics (ASA—Z10.2—1942), prepared by a Committee of the American 
Standards Association, with Society representation, and approved by the 
Association in 1942: 


A 


a 


the area of the cross section, in square feet; 


= the top width of the water surface, in feet; 


_ 1.486 
n 


diameter of the conduit, in feet; 
acceleration produced by gravity; 

the velocity head when flowing full, in feet; 
the velocity head, in feet; 

numerical exponents ; 


A RI, the conveyance of a cross section; 


kp = the conveyance at a full section; 
k, = the conveyance at normal depth; 
Manning’s coefficient of roughness; 


= discharge, in cubic feet per second; 


Q. = discharge when flow is at critical depth; 
Qeayp = capacity of conduit flowing full, with the slope of the 
energy gradient equal to the bottom slope; 


= the hydraulic radius; 


the slope of the energy gradient; 
Spr = the slope of the energy gradient when the discharge is 
flowing through a full section; 
Say = the average value of the slope of the energy gradient at the 
two ends of a reach; 
So = slope of the invert or bottom; 
the average velocity in a section, in feet per second; 
the distance, in feet, measured upstream ; 
a specific reach ; 


= the depth of flow, in feet; 


Yo = the depth of normal flow, in feet; 
¥1,Y¥2 = flow depths at specific sections; 
ye = critical depth; 


= the ratio of the discharge to the capacity of a conduit flowing full; and 
= backwater functions. 


D « 
g => § 
hr = 
im = 
Q 
= 
V => 
z ‘= 
= 
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